On the applicability of an equation of motion method at low-temperatures: 
comments on cond-mat/0309458 and cond-mat/0308413 



o 
o 

(N 
•*-> 
O 

O 



Akira Oguri 

Department of Material Science, Osaka City University, Sumtyoshi-ku, Osaka 558-8585, Japan 

(Dated: February 2, 2008) 

The equation of motion method (EOM) is one of the approximations to calculate transport co- 
efficients of interacting electron systems. The method is known to be useful to examine high- 
temperature properties. However, sometimes a naive application of the EOM fails to capture an 
important physics at low-energy scale, and it happens in recent preprints cond-mat/0309458 and 
cond-mat/0308413 which study a series of quantum dots. These preprints concluded that a unitarity- 
limit transport due to the Kondo resonance, which has been deduced from a Fermi-liquid behavior 
of the self-energy at T = 0, u = [A.O., PRB 63, 115305 (2001)], does not occur. We show that the 
EOM self-energy obtained with a finite cluster has accidentally a singular 1/u) dependence around 
the Fermi energy, and it misleads one to the result incompatible with a Fermi-liquid ground state. 
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PACS numbers: 72.10.-d, 72.10.Bg, 73.40.-c 



I. BACKGROUND 

A Hubbard chain of finite size N connected to two 
nonintcracting leads, which is illustrated in Fig. ^ has 
been studying as a model for a series of quantum dots 
(H 01 [J Of Here N corresponds to the number of 
quantum dots and should be small, while the number 
of the sites in the noninteracting leads at left and right 
are infinite, t.e., Nl — > oo and Nr — > oo. It has been 
shown for the chain of an odd N that the transmission 
probability reaches a unitarity-limit value at T — if the 
system has the electron-hole symmetry (half-filling) to- 
gether with the inversion symmetry (v L — v R ) 4]. This 
conclusion has been deduced from a Fermi-liquid behav- 
ior of the self-energy summarized simply in Eq. J§J . 

Recently, the dc conductance has been calculated with 
an extended version of the equation of motion method 
(EOM) 6], which is essentially the Hubbard I approx- 
imation |3, at zero temperature 0, 0- The results of 
these preprints contradict with the unitarity-limit trans- 
port mentioned above. The purpose of this report is to 
show that the discrepancy is caused just by a naive use 
of the EOM self-energy of obtained with a finite cluster, 
which accidentally has a singularity at low-energies and 
is not compatible with the Fermi-liquid ground state. 

In Sec.|ni we describes briefly the results deduced from 
a Fermi-liquid property Eq. © to make our points clear. 
Then, in Sec. lIIII the properties of the approximate self- 
energy in the EOM are discussed. 
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FIG. 1: Schematic picture of the model: (•) Hubbard chain of 
N sites, and (o) noninteracting leads. The tunneling matrix 
elements vl (vr) connects the cain and the left (right) lead 
to make the energy spectrum of the whole system continuous. 



II. CONDUCTANCE AT ZERO TEMPERATURE 
A. Formulation 

The dc conductance at T = is determined by the 
value of a Green's function G^ 1 (o;) at the Fermi energy 
w = 0H: 
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(i) 



Here GfAuS) is the retarded Green's in the real space 
(1 < ij < N), see Fig. [T] and r a («) = -v 2 Jmg+{u) 
for a — L, R with (g^) being an interface Green's 
function for the semi- infinite lead at left (right). The 
Dyson equation can be written in a N x N matrix form; 



{GH}- 1 = {GV)}- 1 -^ 



(2) 



Here G° = {G^ + } is the noninteracting Green's function 
for U — defined with respect to the total system con- 
nected via tunneling vj, and vr. The inverse matrix of 
the unperturbed Green's function can be expressed as 
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Here 1 is the N x N unit matrix. The self-energy due to 
the Coulomb interaction £ is also defined with respect 
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to the total system of continuous energy spectrum. Note 
that the Green's function Gf N (0) appearing in Eq. 
can be obtained by taking the inverse of the matrix in 
Eq. ©. 



B. Fermi-liquid behavior 



At the temperatures Tr- < T < A gap , we see the insu- 
lating behavior for large TV independent of whether TV is 
even or odd. To observe the Kondo behavior at an ac- 
cessible temperature T, the resonance width Tk should 
be sufficiently large, so that the number of the quantum 
dots TV should be small. 



At T = 0, to — not only the Fermi-liquid but also 
some other "liquid" states of interacting electron systems 
show a behavior, 



ImS(0) = 0, and ReS(O) is finite. 



(6) 



This property of the self-energy, the electron-hole sym- 
metry, and the inversion symmetry are the three keys 
to lock the conductance for odd TV at the unitarity-limit 
value. No other assumptions have been made in the proof 
14 . The perturbation expansion in U of Yamada-Yosida 
|2| preserves the properties of Eq. also for TV > 2 
owing to the presence of the nonintcracting leads. The 
contributions of the hybridization T l and Tr in the un- 
perturbed part G° make it possible to take into account 
the virtual processes in which electrons visit inside the 
nonintcracting leads of continuous spectrum during the 
multi-scattering by U. 

For even TV, the conductance decreases with increas- 
ing TV showing a tendency towards the development of 
a Mott- Hubbard insulator gap A gap Although we 

are interested in the system of small TV, a question arises: 
what can we see in the limit of TV — * oo? Note that the 
results mentioned above hold just at T — 0. To describe 
the thermodynamic limit properly, we must consider the 
physics at finite but sufficiently low temperatures. For 
odd TV, the Kondo resonance of width Tk is situated in 
the Hubbard gap. As TV increases, Tk decreases and fi- 
nally Tk — > at TV — > oo because the developing gap 
disturbs the electrons to form the Kondo singlet. Thus, 
the value of the conductance g N for odd TV (— 2M + I) 
must depend on the order to take the two limit TV — > oo 
and T -> 0: 

• if we take T — > first and then TV — > oo, we see the 
physics at T < Tk, i.e., the Kondo behavior, 



lim 



lim 
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2e 2 /h , 



• if we take TV — > oo first and then T — > 0, we see the 
physics at T > Tk, i.e., Mott-Hubbard behavior, 



lim 



J im 52M+1 

M—*oo 







Note that, for even TV (= 2M), the result does not depend 
on the order to take these two limits, 



,] im 9 2M = 



III. APPLICABILITY OF EOM AT LOW 
TEMPERATURES 

A. Method 

In the EOM calculation of Refs. Q and the self- 
energy has been calculated using a finite cluster of the 
size TVci decoupled from the total system. Typically, the 
cluster introduced for the calculation consists of N inter- 
acting sites, several noninteracting sites in the left lead 
N' L , and that in the right lead N' R ; 



TV- 



N' L 



(7) 



Using the ground state of the decoupled cluster obtained 
with an exact diagonalization, a TYi x TV c i cluster Green's 
function g c \ has been calculated Q, @ . The EOM self- 
energy S c i can be obtained from the Dyson equation for 
the decoupled cluster, 



EdM = {g° cl (u;)} 1 - {?ciMK 



(8) 



Here ^ is the U — Green's function for the decou- 
pled cluster, which describes a discrete energy spectrum. 
Then, an approximate Green's function which takes into 
account the rest of the system can be calculated by sub- 
stituting E c i(w) for the self-energy S(w) defined with re- 
spect to the total system. 

The applicability of this method is determined by the 
separation of the eigenvalues of the cluster Ae. This is 
because E c i has been calculated for the discrete energy 
spectrum, and the correlation effects on the low-lying 
states below this energy separation have not been taken 
into account. Therefore, the application of EOM can be 
justified at high-energies (or high-temperatures) above 
the energy separation Ae. The information about the 
many-body effects on the low-lying states are still miss- 
ing even though one gets a continuous spectrum after 
taking into account the connection to the rest of the sys- 
tem. Therefore, if a fine structure is seen in the obtained 
continuous spectrum, that structure is expected to be 
real only when the energy scale is consistent with the 
resolution Ae. In principle, the energy resolution can 
be improved systematically by increasing the size of the 
cluster successively, and the numerical renormalization 
group (NRG) can be regarded as one of such approaches 
[3, UjJ. Since the cluster is introduced just for the ap- 
proximation and have no physical meanings, the results 
of the observable quantities must not have a strong de- 
pendence on the value of N' L and TVjj. 
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B. The 1/uj singularity of the EOM self-energy 

In Refs. and 0, the calculations have been carried 
out for the cluster with an odd number of sites N c \. At 
half-filling, the ground sate of the cluster is doublet with 
the total spin Stat = 1/2- According to the appendix of 
Ref. 0, the cluster Green's function g c \ is defined using 
the equal- weight sum of S^ ot = 1/2 and —1/2 in Ref. 0, 
while the work in Ref. considers only the cluster state 
with S* ot = 1/2. Furthermore, in Fig. 9 of Ref. 0, these 
two procedures are shown to lead quite different results. 
This means that huge ambiguities exist in the EOM re- 
sults at low energies. In these two procedures, however, 
the equal-weight sum seems to be better because it re- 
stores the rotational symmetry, so that we will use this 
convention in the following discussion of the behavior of 
the EOM self-energy. As typically seen in Fig. 9 of Ref. 
0, the results in Ref. Q have also been affected by the 
singularity, although the position of the pole might shift 
due to the absence of the rotation symmetry of the spin 
in the convention used in Ref. 

The presence of the 1/ui singularity of the EOM self- 
energy for the cluster of an odd size can be deduced along 
the discussion similar to that in the section V of Ref. . 
Using the SU(2) symmetry of the axial charge holding in 
the electron-hole symmetric case (see Appendix 0, the 
cluster Green's function can be expressed as, 



TABLE I: Pole at uj = for the interacting sites of N > 2 



9diii(w) -2t+l 



Ci(l/2,S g + l/2-y\0,S g ;r B ) 



w s Q ( r ') = E i/2,s g ±i/2(r') - E , Sg {r g ) 



x C i (l/2,S fl + l/2;r / |0,5 fl ;r fl ) 

J 1 | } 

\w-<4,(r') w + ws, (»•')/ 

^-y £ £7,(1/2, S g - 1/2; r'|0, S g ; r g ) 

x C 3 (l/2,S g ~l/2-r'\0,S g ;r g ) 

x{ 1 + t^l \, 

(9) 
(10) 



Here S g , J g , and r g are the quantum numbers of 
the ground state, the energy of which is denoted by 
Ej rj .s g { r g)- Note that J g must be zero for U > 0. Fur- 
thermore, for U > 0, the excitation energy must be finite 
W S ( r > m t nc case °f N > 2 as long as the cluster 
size A^ci is finite. Thus, for N >2, the full Green's g c i(tj) 
does not have a pole at ui = independent of whether N c i 
is even or odd. In contrast, the noninteracting Green's 
function g^ { (uj) does has a pole at w = for odd N c i, while 
it does not for even N c \. Therefore the EOM self-energy 
S c i(w), which can be calculated using Eq. ©, must have 
a pole at oj = for odd N c \ because the structure of the 
energy spectrum of g c i(u>) and that of g^{u)) are com- 



Cluster size 


flci(w) 


5&M 


E cl (w) 


Even A^d 


None 


None 


None 


Odd N cl 


None 


Yes 


Yes 



pletely different. On the other hand, for even N c \ the the 
energy spectrum of g c \(u}) and g^^oj) belong to the same 
category having no zero mode, so that £ c i(w) does not 
have a singularity at u> = and shows a behavior con- 
sistent with that of the Fermi liquid Eq. ©. Therefore, 
the presence of the zero-energy pole can be summarized 
as Table [I] (see also Appendix IB1 for N = 1). Note that 
in the case of "Yes" , it does not necessary mean that all 
the matrix elements have the pole. 

Using the EOM self-energy with the Dyson equation 
which connects the cluster to the remaining system, an 
approximate Green's function G^ N can be calculated. Al- 
ternatively, the same G± N can also be calculated by sub- 
stituting the N x N partitioned part of E c i, which corre- 
sponds to the interacting region, into the Dyson equation 
of N x N form Eq. |(SJ|, and then taking the inverse of 
the matrix. Therefore, if one uses the cluster of odd N c \ 
for calculating the value of G^ n (uj) at lu — 0, one obtains 
the result 0; G+ v (0)=0and g N = for odd N because 
of the zero-energy pole of £ c i(oj). In contrast, if one uses 
the cluster of even iV c i in the same situation, the result 
must be different. This also suggests that the EOM is not 
suitable for studying the low-temperature properties. 



C. Interpretation of the zero-energy pole of E c i 

The oscillatory behavior caused by the cluster size A^ c i 
has been seen also in the NRG calculation 

IS 113 
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even N c \, the ground-state is singlet and the self-energy 
can be expanded with respect to u> at the Fermi en- 
ergy w — 0, which ag rees with the local Fermi-liquid 
theory of Nozieres [lj| and the perturbation theory of 
Yamada and Yosida [g. For odd A^i, the ground state 
is doublet, and the self-energy has the zero-energy pole. 
Nevertheless, the energy spectrum of the low-lying states 
shows the correct Fermi-liquid behavior even for odd A^i 
[To| . Furthermore, the residue of the zero-energy pole de- 
creases with increasing iV c i and vanishes in the N c \ — * oo 
limit. Therefore, it is quite dangerous to examine the 
low-temperature properties only by the Green's function 
obtained with a finite cluster of an odd number sites N c \. 

Similar oscillation takes place in an antiferromagnetic 
Heisenberg chain. The spin susceptibility of the chain of 
an odd number sites has a Curie term. Since the weight 
of the Curie term is proportional to l/iV c i and vanishes 
at N c \ — > oo, usually the cluster of an even number sites 
is used for studying the properties in the thermodynamic 
limit, although the cluster of an odd number sites can also 
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be used if the contribution of the Curie term is subtracted 
properly. 

IV. CONCLUSION 

In conclusion, the EOM self-energy obtained for a finite 
cluster of an odd number of sites has the zero-energy pole, 
the residue of which must vanish in the limit of large 
cluster iV c i — > oo. The absence of the unitarity- limit 
transport reported in Refs. |l| and [2j can be explained 
to be caused by a naive use of this singular self-energy. 

APPENDIX A: AXIAL CHARGE 

In the electron-hole symmetric case, the Hamiltonian 
H commute with the total axial charge defined by ^3] 



J 2 


= J 1 , + (J+J_ + J_ J+)/2 , 


(Al) 


Z 


= H 2 (^ViT+^l -1 ) ' 

i 


(A2) 


J+ 


i 


(A3) 


J_ 


= 


(A4) 



Note that the z component corresponds to the charge 
of the U(l) symmetry as Q — 2J Z . The operators J z 
and J± satisfy the commutation relations of the angular 
momentum as the total spin operators S z and S± do. 
Thus, the eigenstates can be classified according to the 
quantum numbers corresponding to J 2 , J z S 2 , and S z ; 

H | J, J z , S, S zV r) = Ej. s {r) | J, J z , S, S z -r) . (A5) 

Eq. (JjJ can be derived using the Wigner-Eckart theorem 
for the spin and for the axial charge, 

(J, J z , S, S z ;r\cl a \J ,J z ,S,S z ;r) 

= (S', S' z ; 1/2, a\S, S z ) (J', J' z ; 1/2, 1/2| J, J z ) 
xCi{J,S;r\J',S'-y) , (A6) 



where the Clebsh-Gordan coefficient appears for the total 
spin (S', S' z \ 1/2, a\S, S z ) and for the total axial charge 
(J', J' z ; 1/2, 1/2 1 J, J z ). Note that the invariant matrix 
element Ci(a'\a) has the following properties against the 
exchange of the arguments a' and a, 

CiG/.S-ijrV-i S;r) 

(A7) 

C i {J,S+~;r'\J-~,S;r) 

(A8) 



APPENDIX B: EOM RESULTS FOR N = 1 

A cluster consisting of a single Anderson impurity at 
the center i = 0, the left lead of N' L sites at i < 0, and 
the right lead of N' R sites at i > can be separated 
into two independent parts by a Unitary transformation 
OLia- = {^Ha + c-ia)/V2 and b ia = {a„ - c- ta )/^/2 for 
i = 1,2,3,... when N' L = N' R (= N^). Only the or- 
bits of die couple with the Anderson impurity to form 
a cluster of the size N' a (= + 1), and the orbits of 
bi a form a noninteracting band of the size N^. Thus, the 
EOM self-energy S c i(w) can be calculated using only the 
cluster of the interacting part, and the properties corre- 
sponding to Table |U can be shown to hold by replacing 
N c \ with N^: £ c i(a/) is zero at lu = and conductance 
reaches the unitarity limit for even N' a , while S c i(w) has 
a pole at ui = and the conductance becomes zero for 
odd N' a . Obviously, the conductance calculated from the 
self-energy having the zero-energy pole contradict with 
the correct unitarity-limit behavior known in the ther- 
modynamic limit N' a — > oo. Note that in Ref. [l|, the 
self-energy for even N' a is used for N = 1, while the self- 
energy having a pole around Fermi energy is used for 
N = 3, 5, 7 (see also Fig. 9 of Ref. @). 
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